An analysis of a regular black hole interior 
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We analyze the thermodynamical properties of the regular static and spherically symmetric black 
hole model presented by Mboyne and Kazanas. Equations for the thermodynamical quantities valid 
for an arbitrary density profile are deduced, and from them we show that the model is thermodynam- 
ically unstable. Evidence is also presented pointing to its dynamical instability. The gravitational 
entropy of this solution based on the Weyl curvature conjecture is calculated, following the recipe 
given by Rudjord, Gr0n and Sigbj0rn, and it is shown to have the expected behaviour for a good 
description of the gravitational entropy. 



I. INTRODUCTION 

The existence of singularities in some solutions of 
, General Relativity is one of the most important unre- 
^ ■ solved issues in our description of the gravitational field. 
, In fact, singularities are an undesirable feature of any 
theory of gravitation: they can be naturallyconsidered 
as a source of lawlessness (see for instance [l|), because 
the space-time description breaks down "there", and 
physical laws presuppose space-time. There are several 
ways in which singularities can be avoided (see [Ij for 
the case of Cosmology), both at the classical and at the 
quantum level. Since there is no widely-accepted theory 
of quantum gravity, the issue of the singularity has been 



addressed from classical considerations in a number of 
ways, one of which is the limiting curvature hypothesis, 
based on the fact that, at a singularity, invariants 
constructed with the Riemann tensor generally diverge. 
Hence, a suitable way to eliminate singularities would 
be to impose a limiting principle on the curvature, asso- 
ciated to a fundamental scale (which may be the Planck 
scale). In practice, this hypothesis (first propposed in 
[sl ) is implemented by imposing that any solution of the 
field equations reduce to a definite nonsingular solution 
(for instance de Sitter space-time) when some of the 
invariants attain their limiting values, in such a way that 
SEC (one of the hypothesis of the singularity theorems 
Hi) is violated. 
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A theory based on this hypothesis has been im- 
plemented by modificating the dynamics of General 
Relativity by the addition of nonlinear terms in Q. 
However, there is another way, suggested by Gliner 
Q, in which de Sitter geometry can be attained for 
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high curvatures, which consists in assuming that the 
microphysics of high-density matter is such that a phase 
transition must occur leading any system under such 
extreme conditions to the de Sitter geometry. Models 
of this type were also considered by Sakharov 9], 
Zeldovich [10], and Bardeen, who suggested that a fluid 
satisfying the condition p = —p could be the final state 
of gravitational collapse [ll|. In fact it was Bardeen who 
in 1968 obtained the first regular solution of Einstein 
field equations having an event horizon, in the presence 
of an electromagnetic field, parametrized by the mass m 
and the charge e. This space-time metric behaves as the 
Schwarzschild metric for large radii and, is de Sitter's 
towards the core of the object. 

Another exact solution of the Einstein field equations 
containing at the core a de Sitter fluid was found by 
Dymnikova [T2j . This solution represents a nonsingular 
black hole, and the stress-energy tensor which is the 
source of the geometry joins the central de Sitter region 
to the standard vacuum at infinity, with an anisotropic 
vacuum layer in between the two. This unsatisfactory 
anisotropic vacuum layer was replaced by Mboyne and 
Kazanas (MK) [Ts'] with a region filled with matter under 
both radial and tangential pressures, described by an 
equation of state (EOS) that relates the radial pressure 
to the density and smoothly reproduces vacuumlike 
behavior near r — tending to a polytrope at larger r, 
low p values Q- 

The main goal of the present work is to study the MK 
model for a regular black hole in which the singularity 
is avoided precisely by the second method mentioned 
above. In particular, a detailed study of the thermo- 
dynamical aspects of the matter which is the source 
of the geometry and also of the gravitational field will 
be given 0. It will be shown that the solution is ther- 
modinamically unstable, and we shall exhibit evidence 
supporting the dynamical instability of the model. The 
issue of the gravitational entropy of this model will be 
also addressed, following the ideas presented in [l6j . 

We shall begin in Sect HI] with a short review of the 
regular black hole solution advanced in [l^. Its thermo- 
dynamical properties will be discussed in Sect lllll and 
the stability will be scrutinized in Sect llVI In SectfVl 
the gravitational entropy of the MK solution will be cal- 
culated. These results will be analyzed in Sect lVIl 



II. A MODEL OF A REGULAR BLACK HOLE 

The model introduced in [l^ is a regular static black 
hole, with a matter source that smoothly goes from a de 
Sitter behaviour near the origin to Schwarszchild's space- 
time outside the object. A spacetime-metric well-adapted 
to examine the properties of this system is 15| : 
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The EOS proposed by Mbonye and Kazanas [iJl is given 
by: 



a 



(a + 1) 



(2) 



where pr is the radial pressure (to be distinguished from 
the tangential pressure, which is needed in these mod- 
els to satisfy the Tolman-Oppenheimer-Volkoff equation, 
p7|). and a — 2.2135 0. Such an equation describes the 
behaviour of matter that changes smoothly from a nor- 
mal behaviour to a core of an "exotic fluid" with an EOS 
that approaches pr ~ — p when r — > 0. The plot of Eq. 
^ is shown in Fig[T] 

In order to solve Einstein's equations, MK chose the 
density profile suggested by Dymnikova |12| : 



p{r) = Pmaxe 

where = SRsTq, Rs is the Schwarszchild radius. 
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Figure 1. Plot of Eq. ([T}, which gives the relation between 
the radial pressure and the density. 



^ Different types of regular black hole solutions were reviewed re- 
cently in [I'M , where regular charged black holes, with the exterior 
described by the Reissner-Nordstrom solution were discussed as 
well. 

^ Other features of the model were studied in llSlI . 



^ This value of a is chosen to yield a sound speed bounded by the 
speed of light [l^ . 
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Figure 2. Pressure as a function of radial coordinate inside 
the regular black hole considered in the text. 



and pmax is of the order of the Planck density. This 
density undergoes a smooth transition from the de Sit- 
ter state at the center to the vacuum state at infinity, 
with an intermediate region of non-inflationary material, 
a situation that was anticipated in [1^. It follows from 
Eq. ([3]) that all the mass is contained inside a sphere of 
radius R within the black hole, which corresponds to the 
surface of the object. 

The plot of the radial pressure as a function of the ra- 
dial coordinate is shown in Fig. [5] It can be seen that 
the pressure follows the equation p = —p at the core of 
the object, and it goes to zero at r/R = 1, the surface of 
the matter region. The radial pressure has another zero 
located at r/R = 0.28, two inflexion points air/R — 0.26 
and r/R = 0.5, and an absolute maximun at r/R — 0.4. 
The exotic matter (for which the pressure decreases with 
decreasing r) occupies the region r/R < 0.40. 

The solution of the Einstein field equations G(j — 
-8nTi/, withTif = diag(p(r),p,(r),pi(r),|7i(r)), for the 
metric given in Eq. ([T|) and the EOS ([2]) takes the form 



r 2 

B{r) = exp / — [m(r') + 47rr'V] 

Jrn ^ 
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2m(r') 



dr' 



where: 



m{r) = An p{r )r dr 
Jo 



(5) 



(6) 



is the ADM mass. 



Although this profile resembles the qualitative plot for the radial 
pressure in a gravastar (see Il7|l ). it must be noted that the MK 
model has a horizon, which is absent in gravastars. 



Equations ([5]) and © describe the geometry of the 
space-time of the regular black hole introduced in p^ . 
Since the EOS and the density as a function of r were 
specified in [is'l, the tangential pressure in terms of the 
given quantities follows from the Einstein field equations. 
Explicitly: 



P± 



Pr + ^Pr 



lip. 



P) 



m{r) + Airr^Pr 
r — 2m(r) 



(7) 



The behaviour of the tangential pressure with r is given 
in Fig. [3l It is important to emphasize at this point 
that the MK model differs from the usual approach (see 
for instance [loj). in the sense that, while in the former 
one the EOS and the explicit form of p{r) were given, 
in the latter two relations between the quantities Pr, p± 
and p are advanced. As shown in the Appendix, the 
choice in the MK model has important consequences in 
the analysis of the perturbations of the field equations. 



III. THERMODYNAMICS OF THE MATTER 

Following [l^, we shall start by assuming that the 
body has reached a static equilibrium configuration (i.e 
inside the black hole, the pull of the gravitational field is 
balanced by the repulsion exerted by the exotic matter). 
It will also be assumed that matter is in thermodynamic 
equilibrium. As shown below, the results obtained 
under these hypotheses turn out to be incorrect, since 
they imply that the system is both dynamically and 
thcrmodynamically unstable. 

The temperature of the matter as a function of the 
radius can be estimated from the laws of standard ther- 
modynamics: 

TdS ^ d{pV) + p dV. (8) 
It follows from this equation, along with Eq. ([2]) and 
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Figure 3. Plot of the tangential pressure as a function of radial 
coordinate, given by Eq. ((7}. 
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Figure 4. Temperature as a function of radial coordinate in- 
side tiie black hole, given by Eq. ((9]). Tsup stands for the 
temperature of the the matter field a,t r = R. 



that: 



sup 



4/3 



,S(r) 



(9) 

where Tgup stands for the temperature of the the matter 
field at r = i?, and: 



S(r) 



ad{p/p„ 



1 + a{p/pn 



(a + l)(p//0max)3 



(10) 

Fig. 13] displays the temperature of the matter as a func- 
tion of the radius. We note that the temperature tends 
to absolute zero close to the core. As for the pressure, 
there is only one maximum at r/R = 0.4. The inflexion 
points of T(r) occur at r/R = 0.31 and at r/R = 0.57. 

An expression for the entropy (up to an additive con- 
stant) as a function of the coordinates can be obtained 
by substituing (O in ©: 



T 



(11) 



The result is plotted in Fig. [Sj The entropy goes to 
zero as r — > in accordance to Nerst's theorem, and 
it has a maximum close to r/R ~ 0.4, in the region of 
highest density of normal matter. The entropy density 
of the matter inside the regular black hole can also be 
calculated, and is given by: 



Sr/2 



(P/Pmax)[l + a{p/PwBx) - (g + l){p/Pu^Bxf] 

(12) 

The plots of the entropy density as a function of p and 
r are shown in Figs. [S] and [T] It can be seen that the 
entropy density diverges at the origin as a consequence of 
the vanishing volume. The entropy density as a function 
of the density has one inflexion point at p/ p^ax = 0.73 
which is located at r/ R = 0.34. 



The radial sound speed as a function of the energy 
density follows from = dp /dp (with c = 1), yielding: 



2{a + l)e-i6'-V«^' 



(13) 



The result is shown in Fig. |8l From Figs. [2] and [H we 
see that the sound speed is zero at the value of r/R for 
which the pressure is maximum. In addition, the sound 
speed takes complex values in the region r/R < 0.4, in 
accordance with the negative slope of the equation of 
state as a function of the density (see Eq. [5] and Fig. [T|) . 
The fact that sound waves do not propagate in the region 
r/i?<0.4isa consequence of the exotic behaviour of the 
fiuid there: a variation in pressure causes an expansion 
rather than a compression. 

To discuss the thermodynamical equilibrium we shall 
need below the Helmholtz free energy, given hy F = U — 
TS = -pV. Explicitly: 



F 



-1 



Fr 



0.2R 



0.018 



(a + 1) 



Pmax 



Pmax 



47r 

Y \R 
(14) 

The plot of the Helmholtz free energy F as a function of 
the radial coordinate is shown in Fig. |9l 

Armed with the functions discussed in this section, we 
shall discuss in the next the issue of thermodynamical 
and dynamical equilibrium. 



IV. EQUILIBRIUM 

IV. 1. Thermodynamical equilibrium 

A condition for a system to be in stable thermodynam- 
ical equilibrium is that, for a given value of entropy and 
volume, the energy must be minimum. This is equivalent 
to impose on the the specific heat at constant volume the 
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Figure 5. Entropy of the matter field as a function of radial 
coordinate inside the black hole. Ssup stands for the entropy 
of the the matter field at r = i?. 
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condition: 



Cv > 0. 



(15) 



The dependence of with the radial coordinate can be 
calculated from: 



(16) 



Using the expressions deduced in the previous sections, 
we obtain: 



V 



For r = R, p is zero, and Eq. pT)) yields: 

^Vsup Pmax 
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sup 



2aT, 



(18) 



sup 



Introducing this expression as a constant of normaliza- 
tion in Eg. pT)) . it follows that: 



Vsup 
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T 
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a-2{a + l) 


{^] 


2' 


T 



(19) 

The plot of the specific heat at constant volume is shown 
in Fig. [TOl The discontinuous behaviour of Cv (with the 
sign change) at r/R — 0.4 is typical of a second-order 
phase transition, suggesting that the system is thermo- 
dynamically unstable. Similarly, we can calculate the 
specific heat at constant pressure, defined by: 



(20) 



p/pmax 

Figure 6. Entropy density of the matter inside the black hole 
as a function of the density. S'sup stands for the entropy of 
the matter field at r = R. 
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Figure 7. Entropy density of the matter inside the black hole 
as a function of radial coordinate. S'sup stands for the entropy 
of the matter field at r = R. 



The result is: 



C'p/Cpsup ^ 



where: 

A = 



/1/2 



a-2(a + l)(^ 



, (21) 



R 



1 



P 



/2 = 1 + 



Pmax 

a - (a + 1) 



a — (a + 1' 
P 



Pn 



(22) 
(23) 



The plot of of the specific heat at constant pressure, 
shown in Fig. [TT] indicates, as in the case of Cy, that 
the system is thcrmodinamically unstable. 

We found that both the specific heat at constant vol- 
ume and constant pressure are not defined for r/R — 0.4, 
where normal matter changes into exotic matter. The 




Figure 8. Radial sound speed as a function of radial coordi- 
nate. 
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Figure 9. Helmholtz free energy as a function of radial coordi- 
nate. Fo.2R stands for the Helmholtz free energy at r/R — 0.2. 
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Figure 10. Specific heat at constant volume as a function of 
radial coordinate in the black hole interior. Cvaup stands for 
the specific heat at constant volume at r = 7?. 
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Figure 11. Specific heat at constant pressure as a function of 
radial coordinate in the black hole interior. Cpsup stands for 
the specific heat at constant pressure at r = ii. 
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Figure 12. Isothermal compressibility as a function of radial 
coordinate, ktr/2 stands for theothermal compressibility at 
r/R= 1/2. 



change of sign of Cv at r/R — 0.4 can be understood in 
terms of Eq. (1161) , the regions where the temperature in- 
creases and decreases, and the density profile given in Eq. 
JH). In Fig. [lUwe can see that Cp = at r/R = 0.57; 
this point coincides with one of the infiexion points of 
the temperature. Cp is positive for r/R < 0.4 and for 
r/R > 0.57. In the region 0.4 < r/R < 0.57, Cp is neg- 
ative. These results show that the specific heats are not 
defined for the same value of r/R where the sound speed 
equals to zero, which reinforces the existence of a region 
of instability for the normal matter field. 

We arrive at the same conclusion by the examination of 
the plot of the the isothermal compressibility kt , defined 
by: 



1 

Kt = 

V 



dV 
dp 



(24) 



as a function of r/R. The equilibrium condition is in this 
case Kt > 0- 

We can summarize the results of this section by assert- 
ing that the discontinuities in the second derivatives of 



the state funcions (along with the continuity of the state 
functions and of their first derivative, as shown in the 
corresponding plots) indicate that the matter inside the 
black hole cannot be in thermodynamical equilibrium. 
This conclusion is reinforced by the plot of the transver- 
sal velocity as a function of the r, defined by: 

2 dp± 

The function v±{r) can be calculated from Eqs. ([2]) and 
([7]) . This function is ploted in Fig. [T31 The plot evidences 
not only the instability discussed above, but also a new 
one, inside the normal matter part of the object. 



IV. 2. Dynamical equilibrium 

To study the dynamical stability of the system, a de- 
tailed analysis of the Sturm-Liouville problem associated 
to the perturbation of the equations of motion for the 
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Figure 13. Tangential sound speed as a function the radial 
coordinate. 



fluid and the metric is mandatory. However, for the rea- 
sons discussed in the Appendix, we shah content here 
with some arguments suggesting that the system is dy- 
namicaUy unstable. First, let us recall that in the region 
of exotic matter there is no propagation of sound waves. 
The equation for the radial sound speed, in the region of 
exotic matter, can be put as follows: 



(25) 



From Eq. (1251) the pressure as a function of the density 
takes the form: 



Ap - 



-vlAp, 



(26) 



where represents the square of the radial sound speed 
and it is always positive. We can see that if the pres- 
sure increases, the density decreases; but if the density 
decreases, the pressure keeps growing. In this process, 
both pressure and density are related in such a way by 
Eq. (pS)) that if the system is perturbed, the fluid never 
stops expanding. We suggest that the huge accumulation 
of energy in this process of continuos expansion might 
lead to divergences that indicate inestability of the sys- 
tem. 



V. ENTROPY OF THE GRAVITATIONAL 
FIELD 

Penrose [20l| suggested that entropy might be assigned 
to the gravitational field itself and he proposed that 
the Weyl curvature tensor can be used to specify it. 
The behaviour of the Weyl tensor follows what is ex- 
pected for the gravitational entropy throughout the his- 
tory of the universe: it is zero in the (homogeneous) 
Friedmann-Robertson- Walker model and it is large in 
Schwarzschild's space-time. 

Rudjord, Gr0n and Sigbj0rn made a recent at- 
tempt to develop a classical description of the gravita- 



tional entropy based on the construction of a scalar de- 
rived from the contraction of the Weyl tensor and the 
Riemann tensor: 



^2 ^ C^^^^C^p^s _ W 



K 



(27) 



This approach is based on matching their description 
of the entropy of a black hole in the event horizont with 
the Hawking-Bekenstein entropy [2l|. In particular, they 
calculated the entropy of Schwarzschild black holes and 
the Schwarzschild-de-Sitter space-time 0. 

Rudjord et al. describe the gravitational entropy of a 
black hole by the surface integral: 



(28) 



where a is the horizon of the black hole and the vector 
field * is: 



* = Per: 



(29) 



where eV is a unitary vector in the radial direction. They 
required Eq. to coincide at the horizon with the 

Hawking-Bekenstein entropy, thus allowing for the cal- 
culation of the constant /cjfl- 

Finally, the entropy density can be determined by 
means of Gauss's divergence theorem, rewriting Eq. 
as a volume integral: 



(30) 



where the absolute value brackets were added to avoid 
negative values of entropy. 

Let us recall that the Weyl tensor is the traceless part 
of the Riemann tensor, and is given by: 

2 

H 7^{9a[iR5]l3 - 9l3[iR5]a) + (31) 

2 

+ (n-l)(n- 2)^^-1^^^!^' 

where Rap-yS is the Riemann tensor, Rap is the Ricci ten- 
sor, R is the Ricci scalar, [ ] refers to the antisymmetric 
part, and n is the number of dimensions of space-time. 

The absence of structure in space-time corresponds to 
a null Weyl conformal curvature {W = C°-^'"^Cabcd = 0). 
Hence, the Weyl tensor contains the information of the 
gravitational field in absence of matter and other non- 
gravitational fields. 

The Weyl scalar for the space-time metric given in Eq. 
© is: 



2 

max V 



PiX 

ttR^Pu 



P3X +P2X + 

.x)r^R^], (32) 



See reference 12211 for other space-times 
^ To simplify the notation, we will set the constant ks equal to 1 
since it plays in what follows the role of a scale factor. 
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Figure 14. Weyl scalar as a function of radial coordinate. 
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Figure 15. Kretschmann scalar as a function of radial coordi- 
nate. 



where x — e 
?6 



'^"r^ and the coefficients pi are: 



- 192r6a^i?3pmax - 2R^r^aTTp^^, 
P3 = 2i?^r^a7rp,„ax + 144r^a7ri?3/9max, 
P4 = (-1152 - 1152a)r^ - 487ri?Vmax(-8a - 8 + a^)r^ + 

+ 2TTp^,^R%a + iy, 
P5 = ~3367rp„,axi?^(a + 1^ - 27rp„,axi?^a + 1^, 
P6 = 967rp„iaxai?^(a + l)r^, 

PS = -487ri?3p^^^^(a + 1)2^6. (33) 

If we let r — ^ 0, the Weyl scalar goes to zero. Since the 
regular black hole space-time has a de Sitter geometry in 
the origin, this is the expected limit. 

The plot of the Weyl scalar is shown in Fig. [141 We see 
that for large values of r the Weyl scalar tends asymp- 
totically to zero. In the matter region {r/R < 1) it has 
one absolute maximum at r/R = 0.5 and one relative 
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Figure 16. Gravitational entropy density as a function of ra- 
dial coordinate. 



maximum at r/R = 0.26. Between these points, 14^ = 
at r/R — 0.34. If we analyse the equation of state as a 
function of the density (Equation (2), we find one inflex- 
ion point at r/R ~ 0.5. The absolute maximum of the 
Weyl scalar seems to be related to the transition point 
between the two regions of matter. 

The value oi r/R for which the Weyl scalar has one 
relative maximum is close to the point where the pres- 
sure is zero. This suggests a relation between the re- 
gion where matter has a negative pressure and the be- 
haviour of the Weyl scalar. Notice that the inflexion 
point of the entropy density of the matter coincides with 
the value of r/R for which the Weyl scalar is zero, that 
is at r/R — 0.34. At this point the gravitational field 
changes from attractive to repulsive. 

We also calculate the Kretschmann scalar. The plot is 
shown in Fig. [15] Again, we see that outside the matter 
region the Kretschmann scalar tends asymptotically to 
zero. If we let r ^ 0, it goes to: 



K = 



5127r2pn 



(34) 



as expected for the de Sitter space-time. The 
Kretschmann scalar is always positive and it has one ab- 
solute maximum at r/R = 0.27 and one relative maxi- 
mum dXr/R — 0.47. 

As a cross-check, we see from the plots that the calcu- 
lated Weyl and Kretschmann scalars satisfy the relation: 



(35) 



that holds in every static spherically symmetric space- 
time [l^. The gravitational entropy density of the MK 
model can be computed following following Rudjord et. 
al.'s proposal. The plot of the entropy density is shown 
in Fig. [TH It is seen that s = at r/R = 0.27, which 
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is close to the transition point from positive to negative 
pressures. The absolute maximum is at r / R — 0.37 and 
the two relative maxima are at r/R — 0.16 and at r/i? = 
0.65. For large values of r/R the gravitational entropy 
density tends asymptotically to zero and in the core of 
the object it is also zero. This last result is in accordance 
with a correct classical description of the entropy of the 
gravitational field. 



VI. CONCLUSIONS 

We have shown that the thermodynamical quantities 
describing the matter that is the source of the regular 
black hole model proposed by Mbonye and Kazanas in- 
dicate that the region of exotic matter is unstable. Some 
evidence has been presented that points to the dynam- 
ical instability of the model. The evidence for this sec- 
ond instability is supported by the plots of Vr and 
In particular, the latter indicates the existence of a sec- 
ond type of dynamical instability, in the normal matter 
part of the object. These findings should be confirmed 
by a perturbation analysis (which will be presented else- 
where), based on the equations given in the Appendix. 

We have also calculated the Weyl and Krestchmann 
scalars for the regular MK black hole, and a possible 
classical estimator of the gravitational entropy proposed 
by Rudjord, Gr0n and Sigbj0rn based on the Weyl cur- 
vature conjecture. It was shown that close to the core 
of the object, the entropy density tends to zero as well 
as for large values of the radius. Hence, this classical 
estimator gives a good description of the entropy of the 
gravitational field in de Sitter and Schwarzschild limits. 



APPENDIX: DYNAMICAL STABILITY OF THE 
INTERIOR REGION 



In order to study the stablility of a system with 
tangential pressures against radial perturbations, both 
Einstein field equations and the equation of state must 
be perturbed following the standard procedure (see 
for instance [ISl)- The goal is to obtain a differential 
equation for the radial dependence of the quantity ^ (r, t) , 
which represents a small radial displacement of the fluid 
respect to its position of equilibrium at time t, that is, 
rir,t) = ro +^(r,t). 

In what follows, the unperturbed metric coefficients 
A(r), B{r), and the unpertubed thermodynamics vari- 
ables will be denoted by a zero subindex. The corre- 
sponding perturbations are denoted by SA(r, t), 5B{r, t), 
Spr{r,t), 6p±{r,t) and Sp{r,t). A long but straightfor- 
ward calculation (which parallels that presented in [23|) 
shows that the equations for the perturbations are given 



by: 

Snr'^Sp 
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SB\ r Aq d^^ 
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Notice that there are six unknowns in this system, a 
consequence of the choices made by MK (one EOS and 
the explicit form of p{r) to determine the system). 

Equation can be integrated respect to time Q, 
yielding: 



5A 

^0 



-87rr(po-|-rf)AoC. 



By isolating 5B from Eq. ([57]) wc get: 

d fSB\ „ ^ ^ 5A / 1 dBa 

The final set of equations takes the form: 

d fr6A' 
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Standard manipulations using the expression ^ (r, t) = 
a{r) exp iut in these equations lead to the following sec- 
ond order ordinary inhomogeneous differential equation 
for the function (T(r): 

A(r)a"(r) + a(r)a'(r) + [B^r) + ujD^r)] a{r) = n(r), 

(43) 

where we have also used Sp±{r, t) = n(r) exp iujt. 

This equation defines an inhomogeneous Sturm- 
Liouville (SL) problem, differently from the common 



We choose the constant of integration so that &A(r, t) = if 
&r = 



10 



case which yields a homogeneous one, the inhomogeneity 
being a direct consequence of the way the MK model 
solution was obtained, as mentioned before. 

Equation p3)) is to be solve for r G [0,1], using 
the boundary conditions cr(r) — 0, (7(1) — 0, and 
[6p{l,t) + p'^{l,t)£_{l,t)] = [23|. The coefficients A^{r), 
B^{r), Ci,{r), and D^,{r) are given by: 



(r) = Bi (r) + (r) + (r) + E4 (r) , (45) 
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where: 



Figure 17. Coefficient A* as a function of radial coordinate. 
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For a given n(r), Eq. P5|) might be solved using the 
Green's function method (see for instance |23]). However, 



as we shall show next, the SL problem defined by Eq. p3)) 
has a singular point. Consider the equation [25i] : 

- ipy')' + qy^ Awy, (55) 

which corresponds to the homogeneus SL problem de- 
fined on on J = (a, &), with —oo < a < b < oo. The 
eigenvalue A is such that A G C, and p,q,w and y are 
functions of x. It is also assumed that: 



1 

7 

P 



q,w e iioc {J,C) , 



(56) 



where Lioc (J, C) denotes the linear manifold of functions 
y satisfying y £ L([a,/3],C) for all compact intervals 
[a, f3] C J. L ( J, C) denotes the linear manifold of com- 
plex valued Lebesgue measurable functions y defined on 
J for which: \yit)\ dt = Jj \y{t)\dt = Jj \y{t)\ < 
Here J is any interval of the real line, open, closed, half 
open, bounded or unbounded. Following [25|, we have 
the next definitions: 



• The (finite or infinite) endpoint a is regular if 
1 



P 



-,q,w e L{{a,d), 



(57) 



holds for some (hence any) d ^ J. 



• An endpoint is called singular if it is not regular, 

with similar definitions at r = 6. To analyse whether 
a = is a regular or singular endpoint for our problem. 



defined by Eq. (^5]) . we show in Fig. [T7] the plot of 
the coefficient A^, as a function of radial coordinate. We 
can see that A^{0) = A^{OA) = 0, A^{r) < for r G 
(0,0.4), and A^r) > for r G (0.4,1). Hence, ^ 
L((0,c?),C) for any d G (0,1). Therefore, the endpoint 
a = is singular. If an endpoint is singular it could be 
a limit point or a limit circle. According to [25^, "there 
is hardly any literature on the LP/LC (limit point-limit 
circle) dichotomy when all three coefficients are present 
in the SL equation". In particular, "there seems to be no 
literature on LP/LC criteria when p changes sign". Due 
to these complications, we shall attack this problem in a 
future publication. 
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